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Wc find that generic entanglement is physical, in the sense that it can be generated in polynomial 
time from two-qubit gates picked at random. We prove as the main result that such a process 
generates the average entanglement of the uniform (unitarily invariant) measure in at most 0(N 3 ) 
steps for N qubits. This is despite an exponentially growing number of such gates being necessary 
for generating that measure fully on the state space. Numerics furthermore show a variation cut-off 
allowing one to associate a specific time with the achievement of the uniform measure entanglement 
distribution. Various extensions of this work are discussed. The results are relevant to entanglement 
theory and to protocols that assume generic entanglement can be achieved efficiently. 

PACS numbers: 03.67.Mn, 05.70.-a 



Introduction — Entanglement has traditionally been 
viewed as a fundamental tool for studies of the founda- 
tions of quantum mechanics More recently, the view- 
point of using entanglement as a resource has also gained 
prominence; see Q for a recent review. While a great deal 
of insight into the structure of two-particle entanglement 
has been gained, it has become equally clear that the 
complexity and diversity of multi-particle entanglement 
grows exponentially with the number of particles. It is 
thus difficult to imagine a structurally simple theory that 
characterizes and quantifies all details of multi-particle 
entangled states. On the other hand one may expect 
that large numbers of particles admit a notion of typical 
entanglement properties for which a structurally simple 
theory may be developed. This intuition gives hope that 
significant progress can be made by restricting attention 
to entanglement properties that are typical (generic) rel- 
ative to the uniform (unitarily invariant) measure, the 
unbiased distribution of pure states. In this setting it was 
demonstrated that typically pure states of large numbers 
of spins exhibit maximal bi-partite 0, 0, [H, Q and multi- 
partite entanglement This suggests that the explo- 
ration of the entanglement properties of generic states is 
a promising approach. 

But a big question mark exists as to whether statements 
about generic states relative to the uniform measure are 
physically relevant. This is because the generation of 
a typical unitary requires a sequence of 2-qubit unitaries 
whose length grows exponentially in the number of qubits 
0, even if one allows for a finite fixed fidelity. Thus 
achieving the uniform distribution to a fixed accuracy 
requires sequences of random 2-qubit unitaries that grow 
exponentially with the size of the system, and quickly 
becomes unphysical, see e.g. d, H, llO, lU|. One could 



then argue that the entanglement properties of generic 
quantum states are mathematically sound and interest- 
ing but physically irrelevant, as a system undergoing a 
randomisation of its state through two-party interactions 
would only get close to the uniform measure in an unfea- 
sibly long time. On the other hand, entanglement prop- 
erties represent a restricted class of physical properties of 
a quantum state. Accordingly the faithful reproduction 
of generic entanglement properties may be possible with 
far fewer physical resources, i.e. 2-qubit gates, than those 
required for the generation of the expectation value for 
an arbitrary observable. 

It is thus crucial to explore whether generic entangle- 
ment properties can be obtained efficiently, i.e. poly- 
nomially in the number of qubits, using only one- and 
two-qubit gates. The present work answers this question 
positively (c.f. Figure I). 
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FIG. 1: Typical numerical simulation using the random cir- 
cuit. The entanglement average of the uniform measure is 
reached to an accuracy e in n steps. We prove that it suffices 
with n = 0(N 3 ) to achieve a fixed arbitrary e accuracy when 
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ploration of generic entanglement towards a structurally 
simple entanglement theory and have direct practical rel- 
evance since certain quantum information processing pro- 
tocols such as [1, [lU, [IH, [HI assume that generic entan- 
glement can be generated efficiently. 

The presentation proceeds as follows. We firstly define 
the key process that is used throughout this work: ran- 
dom two qubit interactions, modeled as random circuits 
on a quantum computer. Then we prove that the generic 
entanglement average as well as the purity of a subsystem 
are achieved efficiently and that the so generated states 
are typically very close to maximally entangled. This is 
followed by numerical evidence that the achievement of 
generic entanglement can be associated with a specific 
time, the variation cut-off, for large systems. We finish 
with a discussion and conclusion. 

The setting - We consider a set of iV-qubits split 
into two subsets A (with N A qubits) and B (with Nb 
qubits). Let | tpo) be a initial state in AB and consider a 
random circuit C n consisting of n randomly chosen two- 
qubit quantum gates. Define \ip n ) — C n \ipo) and the 
reduced density matrix p An — Trg(\ ip n ) (^n I) of sys- 
tem A. Then the entanglement in the state is given by 
E(i> n ) = S(p At „) and its purity by Tr(p 2 A n ). 

Definition of the random circuit - The random 
circuit C n is a product W n . . . W\ of two-qubit gates 
where each W\ is independently chosen in the following 
way: A pair of distinct integers c 7^ t is chosen uniformly 
at random from {f , . . . , N}. Next, single-qubit unitaries 
U[c] and V[t] acting on qubit c and t respectively are 
drawn independently from the uniform measure on U (2). 
Then W = CNOT[c,t]U[c]V[t] where CNOT[c,t]js the 
controlled-NOT gate with control c and target t |2ll ]. 

Asymptotics of random circuit — The circuit, act- 
ing on N qubits, will asymptotically induce the uniform 
measure on states (see e.g. [§[). In the above setting 
for states distributed according to the uniform measure 
the average bipartite entanglement can be found exactly 
@ and is bounded from below such that E[E(i() n )] > 
N A - T ^2~ t where N B -N A = t > 0. Likewise we find 
that the average purity of the subsystem A is given by 
(2^+2^)7(2^ + 1) consistent with [13]. Furthermore, 
the distributions for entanglement and purity concentrate 
around their average with increasing N [!, 0, H, [1] . Thus 
one is overwhelmingly likely to find near-maximal entan- 
glement for large systems. 

Main Theorem — We will now be concerned with the 
approach to the asymptotic regime. For the above set- 
ting we prove that, independently of the initial state IV'o), 
convergence of the expected entanglement to its asymp- 
totic value to an arbitrary fixed accuracy e is achieved 
after a number of random two-qubit gates that is poly- 
nomial in the number of qubits. More precisely we find: 

Theorem 1 Suppose that Nb — N A = t > and that 
some arbitrary e <E (0, 1) is given. Then for a number n 
of gates in C n satisfying 

n > 9N{N- l)[(31n2)A + lne" 1 ]/4, 



we have E [E(ip n )] > N A — (2~* + e)/ In 2. (1) 
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Eq. ([2]), follows from cq. (JTJ) employing y/2S(cr\\p) > 
tr\a — p\x, where S is the relative entropy, and \tr\a — 
p\x > 1 — tr^J ^fpa^fp as well as Uhlmann's TheoremQ. 
To prove eq. |TJ) we prove a Lemma that considers the 
quantity EfTr^J]. 

Lemma 1 For arbitrary TV, N A , Nb and all n we have 
E [Tr{p\ n )) - - 



<A N e 



2 N + 1 

To see that Lemma [T] implies Theorem [T] note first that 
E(ij) n ) = S{p An ) > -\og 2 Tr(p 2 An ). By convexity 
we then find E [- log 2 Tr(p 2 A n )] > - log 2 (E [Tr{p% n )] ) 
and a direct computation using ln(l + x) < x for x > 
completes the argument. 

Proof of Lemma 1 - We proceed to outline the 
proof of Lemma Q] below, omitting some tedious but 
straightforward calculations to improve clarity. Wc 
begin with a useful representation of quantum states 
in terms of Pauli-operators. Indeed, |VVt)(VVi| = 



E 

2 -N/2 Tr( 



P6{0,: 



•£»(p)2" 
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a Vi [i] , where each (p) = 
(M> n |) and a Pi [i] is a Pauli oper- 



ator acting on qubit i [21| . Then for the reduced density 
operator p An = Tr B {\ i> n ){ip n I) we find 



E [Tr(pXn)] 



■>N, 



(3) 



The main purpose will now be to analyze the evolution of 
the expected values of the squared coefficients, E [£^(p)l 
Evolution of the coefficients - The key idea of 
the proof relies on the observation that the E [£^(p)] 
form a probability distribution on {0,x,y, z} N for all n 
and that these probabilities evolve as a Markov chain 
with transition matrix P which takes q distributed 
according to (E [£n(<z)] )<? m onc s t e P to p distributed 
according to (E [£n+i(p)] )p- To determine P we con- 
sider the action of a random unitary W n at time n 
that acts on qubits c, t in state \ip n }- This results in 

E [&+i(p)\ipn,c,t] = 1/16 E £n{q)£ n W)x 

q,q'e{0,x,y,z} N ^i^{c,t},q i ^q' i 

xE \Tr\Uo-i- [c}WaPc [c}}Tr[Ucr< [c]W(jP- [c 

xE Tr[V<7 qt [t}V^aP t [t}}Tr[Vait[t}V''aP t [t}} 

where the (p c ,pt) are uniquely determined by 
CNOT[c, t]o-P° [c]cr pt [t]CNOT[c, t]=±aP° [c]a^ [t]. Direct 
calculation with the uniform measure on U(2) shows that 
the products of expectations in the sum vanish unless 
q = q' . Then with the Kronecker symbol Sij we find 

E[e n+1 {p)\^n,c,t\ = e p^ t] {q,p) n s q ^e n (q) 

qe{0,x,y,z} N i£{c,t} 
where P^(q,p) = l if Pc =p t =q c =q t =0- P^ 4 ) (q, p)=l/3 
if Pc = q c = and p t ,q t ^ or if p t = q t = and 
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PcQc^O; and P' c,t '(q,p) — 1/9 otherwise. Averaging 
P^(q,p) n^{ c . t } $qi., P i over the N(N-l) choices of c, £ 
produces the entry P(q,p) of the transition matrix of 
the desired Markov chain for (E 

Simplifying the Markov chain - Our aim is the 
evaluation of eq. ([3]) and it turns out that this can be 
done via a simplified Markov Chain. Consider {q(n) = 
(q(n)iq(n)2 ■ ■ ■ 9(^)jv)}n>i as an n-step evolution of our 
Markov chain P. Then the sets S(n) = {i € {1, . . . , N} : 
q{n)i jtz 0}, identifying the nonzero elements of q(n) also 
form a Markov Chain. Using S(n) in eq. we find 



E [Tr{p% n )] =2 w "P(£(n)c A). 



(4) 



Thus we need only to consider the chain {S(n)} n . 

Convergence rate of the Markov chain. - As 

it turns out, our chain is almost ergodic: removing 
the isolated state S(n) = 0, we obtain an ergodic 
chain on = 2* 1 -- A '>\{0}. Since P(S(n) = 0) = 
P(5(0) = 0) = 2~ N , determining the convergence rate 
to the equilibrium of S(n) on Cl given by M(S) — 
3ls|/( 4 JV _ ^ s g Q ig sufficient for 

our purposes. Let 

Q = {Q{S, S'))s,s'en be the transition matrix of the re- 
stricted S(n) chain. It has largest eigenvalue 1 whose 
eigenvector determines the steady state solution M. . The 
difference to the second largest eigenvalue, the spectral 
gap Aq , bounds the convergence rate to the steady state: 
for any initial distribution vector v, the component of 
Q n v orthogonal to M. shrinks exponentially fast with 
Xqu. A quantitative result is provided in Chapter 2 
of [24[ (see Corollary 2.15): since our Q is a reversible 
chain with Q(S,S) > 1/2 for all S € Cl, we obtain 
|P (5(n)c Ayj2^ ScA M(S)\ < e- A «"/\Amn T Al(T) < 

2 JV e" A « n . We have Ee^scA M ( s ) = - V/i^ ~ 
1). Putting back the isolated state into the calcu- 
lations, applying (J4j) and noting that 2 B < 2 N yields 



2" A +2' 



2 N + 1 



< 4 N e 
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All that remains 



E [Tr{p% n )] 

is to show that Aq > 4/9N(N— 1). We use a well-known 
variational principle for Aq [251 ] : 

_ . Z s . s . € nM(S)Q(S, S')(f(S) - f(S')) 2 
9 m ET,T>enM(T)M(Ti)(f(T)-f(T>)r> U 

where the inf is taken over non-constant / : Cl — > R. 
This is an application of Raleigh's principle to the sec- 
ond smallest eigenvalue of I — Q, which is precisely Aq. 
Eq. ([5]) implies that if R is the transition matrix of a 
Markov chain on CI with same stationary distribution A4 
and aR(S,S') < Q(S,S') for all S,S' € Cl, then the 
gap \r of R satisfies Aq > aXji. This allows us to es- 
timate Aq by comparison with a simpler chain [23| . In- 
deed, our R will be transition matrix of chain {B(n)} on 
Cl defined as follows: Assume B(n) = B and choose a 
1 < j < uniformly at random. If j e £> and \B\ > 2, 
set B(n+l) = B\{j} with probability 1/3 and B n+ i = B 
with probability 2/3. If j € B and \B\ = 1, do noth- 
ing. If j £ B, set B(n + 1) = B U {ff. This is a 
biased random walk on the hypercube 2^ 1, " ,JV '' where 



transitions to state are suppressed. A coupling argu- 
ment following e.g. Chapter 4 of [12] shows that R has 
a spectral gap > 1/3N. Moreover, one can check that 
aR(S,S') < Q(S,S') with a = 4/3 (N - 1). It follows 
that Aq > aXn > 4/9A r (A r — 1), as desired, and the 
proof is finished. Numerics indicate convergence in ap- 
proximately NlogN steps, so our bound is not tight. 

Observe cut-off - Many Markov chains exhibit the 
so called " cut-off effect" [20] . The cut-off refers to an 
abrupt approach to the stationary distribution occurring 
at a certain number of steps taken in the chain. Say 
we have a Markov chain defined by its transition matrix 
P, and that it converges to a stationary distribution 7r. 
Initially the total variation distance TV —\\ P — ir \\ = 
sup | P(E) — n(E) | between the corresponding probabil- 
ity distributions is given by TV — 1. After k steps this 
distance is given by TV(k) =|| P k — n \\. A cut-off occurs, 
basically, if TV(k) ~ 1 for k = 0, 1,2, ...a and therafter 
falls quickly such that after a few steps TV(k) ~ 0. As 
we increase the size of the state space, the ratio of the 
number of steps during which the abrupt approach takes 
place and a should vanish asymptotically. Then we can 
say that the randomisation occurs at a steps. Rigorously, 
this may be stated as follows 20]. Let P n , n n be Markov 
Chains on sets Xn- Let a n , b n be functions tending to in- 
finity, with b n /a n tending to zero. Say the chains satisfy 
an a n , o n cutoff if for some starting states x n and all fixed 
real 9 with k n = \_a n + 9b n \, then || P%" — tt„ j| — > c{9) 
with c(9) a function tending to zero for 9 tending to in- 
finity and to 1 for 9 tending to minus infinity. Here we 
observe this behaviour in the entanglement distribution, 
a functional of the Markov chain on unitaries given by the 
random circuit, and we accordingly term this a cut-off. 
Numerical Observation- Numerical simulations indi- 
cate a cut-off effect in the entanglement probability dis- 
tribution under the random circuit on \ fy® N may be ob- 
served both for single qubit gates drawn from the uniform 
measure on U(2) and for stabilizer gates; see Figure^ 

The simulations using stabilizer gates allow us to con- 
sider far larger systems sizes. Here we choose the single 
qubit gates U and V from the set {a x , a v ,a z , S, H} [Hj 
with equal probability. It should be noted that the proof 
of Lemma 1 still holds (see [27| for details) and the en- 
tanglement behaviour will remain similar [28|. The re- 
striction allows us to use the efficient stabilizer formalism 
[TtJ and the tools developed in [l8| which in turn allows 
for an efficient evaluation of state properties. 

Extensions of the present result - Similar meth- 
ods can be used to address the mixed state setting 
through tracing out part of the system on which the 
random circuit is applied. Multipartite entanglement 
measures based on average purities [2j can be considered 
with the results established here. We also anticipate that 
one can use similar techniques to obtain rigorous state- 
ments about the convergence rates of finite temperature 
Markov Process Quantum Monte Carlo simulations. 

Our results may be applied to the protocols for super- 
dense coding of quantum states presented in 0, [l2|, [HJ 
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FIG. 2: Observe a variation Cut-Off of the entanglement 
probability distribution compared with that of the uniform 
measure as determined numerically. The state space has been 
discretised by rounding off entanglement values to the nearest 
integer. We observe that TV — 1 for a while and then falls. 
Finally there is a stage where TV ~ 0. The effect becomes 
more pronounced with increasing N. The results for N > 8 
are done using the stabilizer random circuit. 

to replace the inefficient process of creating random uni- 
taries distributed according to the Haar measure by our 
efficient random circuits. After that replacement, Theo- 
rem I may be ap plied directly to verify that the main 
Lemma 1 of [1J| still holds. It is an open ques tion 
whether the performance of the protocols in pj 1 12L Il3j 



is adversely affected by this substitution. This cannot 
be decided on the basis of Theorem 1 alone but we ex- 
pect that similar techniques as described here and in [27| 
will be able to decide this. The results of this work as 
well as the above extensions will be presented in detail 
in forthcoming publications. 

Conclusion — In this work we have proved that 
the average entanglement over the unitarily invariant 
measure is reached in a time that is polynomial in 
the size of the system by a quantum random process 
that is restricted to random two-qubit interactions. We 
also provided numerical evidence that for large systems 
the entanglement distribution of the uniform measure is 
achieved at a specific point in time, the variation cut-off. 
Our results demonstrate that the entanglement proper- 
ties of generic entanglement are physical in the sense that 
they can be generated efficiently from random sequences 
of two-qubit gates. We have described extensions, includ- 
ing how this knowledge can be applied to render certain 
protocols efficient. 
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